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Abstract: The measurement of the low transverse momentum region of vector boson pro- 
duction in Dreh-Yan processes has long been invaluable to testing our knowledge of QCD 
dynamics both beyond fixed-order in perturbation theory as well as in the non-perturbative 
region. Recently the D0 collaboration have introduced novel variables which lead to im- 
proved measurements compared to the case of the standard Qt variable. To complement 
this improvement on the experimental side, we develop here a complete phenomenological 
study dedicated in particular to the new (j)* variable. We compare our study, which con- 
tains the state-of-the-art next-to-next-to-leading resummation of large logarithms and a 
smooth matching to the full next-to-leading order (O (a^)) result, to the experimental data 
and find excellent agreement over essentially the entire range of (j)* , even without direct 
inclusion of non-perturbative effects. We comment on our findings and on the potential for 
future studies to constrain non-perturbative behaviour. 
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1. Introduction 

The transverse momentum (Qt) distribution of lepton pairs, or equivalently vector bosons, 
produced via the Drell-Yan process [1] has been a classic observable in the realm of phe- 
nomenological studies of QCD at hadron colliders, both in and beyond the perturbative 
domain with pioneering studies commencing over three decades ago [2-7]. In spite of be- 
ing well-studied, the Drell-Yan Qt variable remains of continuing importance both on the 
QCD theory side as well as for high precision Standard Model phenomenology, such as W 
mass determination, and hence with important implications for Higgs studies as well as 
potential studies involving new particles that may be discovered for instance at the LHC. 

As far as QCD theory is concerned, while the high Qt tail of the distribution ought 
to be described within fixed-order perturbative methods, the low Qt region is enriched by 
the presence of large logarithms which require resummation in order to yield a meaningful 
result. A successful description of the entire Qt distribution hence requires resummation 
(at least to next-to-leading (NLL) accuracy) accompanied by matching to, ideally, NLO 
fixed-order results. In fact the state of the art for the Qt variable as far as resummation 
is concerned is up to the NNLL level [5,8,9] and fixed-order codes such as MCFM [11], 
DYNNLO [12] and FEWZ [13] offer the required fuh NLO resuhs for the differential dis- 
tribution. Thus, on the theoretical side, the level of accuracy one can achieve in the study 
of the Qt variable is paralleled by few other observables in QCD. Successful confrontation 
of such calculations with data indicates control over QCD dynamics at the high precision 
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level and paves the way for extending such studies to related observables at hadron col- 
liders, such as the Higgs Qt distribution [14]. There are also relatively new calculational 
approaches such as those of soft-collinear effective theory where the computation of the Qt 
spectrum has been carried out [9,10,15] and yields comparable results to traditional resum- 
mation approaches as well as yielding the full calculation of the next-to-next-to-leading 
logarithms, a piece of which had not previously been computed till date [9]. 

Armed with such accurate perturbative predictions one should be ideally placed to 
explore the role of non-perturbative (NP) effects and either extract them from the data or at 
the very least set limits on their size. One may visualise the NP behaviour as, for instance, 
an intrinsic Fermi motion of partons within the proton which leads to a Gaussian smearing 
of the perturbative Qt spectrum, with comparisons to data offering an opportunity to 
constrain the parameters of the Gaussian. A commonly used parametrisation of non- 
perturbative effects is the Brock-Landry-Nadolsky-Yuan (BLNY) form factor [16] which 
was obtained by comparing the resummation formalism of Ref. [7], as encoded in the 
RESBOS generator [17, 18], to Tevatron Run-I data. 

Alternatively, one can take the perturbative calculation as a means of directly ap- 
proaching NP behaviour and study the ambiguities inherent in the perturbative approach 
(due to the divergence of the perturbative coupling). For examples of such studies we 
point the reader to Ref. [19] for the Qt case and the detailed exposition in Ref. [20] for the 
related case of energy-energy correlation in e'^e~ annihilation as well as references therein. 

It is perhaps worth noting here that various approaches to Qt resummation do not 
concur on the size of non-perturbative effects required to describe experimental data. As an 
instance of this we note that in Ref. [21] an essentially perturbative approach was seen to 
give an adequate description of Tevatron Run-II data except at the very lowest Qt values 
where one may expect very strong NP effects. As a general remark, in the above context, 
before reaching firm conclusions on the size of NP effects (which one should in any case 
expect to vary from observable to observable and depend on the process) one should be sure 
about the robustness of the perturbative result and the uncertainty associated to it, which 
can be ascertained by varying the various perturbative scales in the problem. We shall 
comment in more detail on the role of these scale variations later in this article. Lastly, 
we should also point out that given the sensitivity of the Qt spectrum to QCD dynamics 
beyond fixed-order, it is an ideal variable for the testing of Monte Carlo event generators 
and for tuning the parameters therein as well as for testing new models of non-perturbative 
behaviour [22] encoded in QCD event generators. 

Recently, the D0 collaboration have introduced two new variables o/r and (f)* which 
have been shown to offer advantages over the standard Qt variable in that they can be 
measured more accurately and hence offer the potential to push the theoretical studies 
further, yielding valuable new information on both perturbative and non-perturbative ef- 
fects [23-25]. As a concrete example of the utility and discriminating power of the new 
variables one can consider the issue of the small-x broadening of non-perturbative effects 
suggested initially in comparisons of the RESBOS generator to semi-inclusive DIS data [26]. 
Since such small-x effects could have a serious impact on the vector boson and Higgs re- 
summed Qt spectra, especially in forward (high rapidity) regions and, in particular, at 
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the LHC, it is clearly of importance to ascertain their presence or otherwise. Based on 
comparisons to their new data on cp* the D0 collaboration were able to demonstrate that 
such effects were in fact disfavoured, which had not been previously possible due to errors 
on the Qt spectrum even with Tevatron Run-II data [25]. 

In previous papers [27, 28] we have provided the details of a resummed treatment of 
the new variables and discussed their relationship to Qt and to each other. We have 
computed the resummation to NNLL accuracy and carried out the matching to fixed-order 
NLO results from MCFM. Since the variables ax/M and (p* are essentially identical at low 
Qt [24,25,28] we shall focus here on the (f)* case which is also the variable favoured by the 
D0 collaboration in terms of measurement [25] . In the present paper we extend our initial 
theoretical studies to the phenomenological level by deriving our matched resummed results 
differentially in the vector boson rapidity and with precisely the same cuts as adopted for 
the D0 measurements [25]. We study the role of various scales in the problem including 
renormalisation and factorisation scale variation as well as resummation scale uncertainties 
(rescaling the argument of the logarithms we are resumming) and derive uncertainty bands 
for the perturbative result. We then compare our results to the experimental data for both 
electrons and muons in various rapidity bins. 

We organise the present paper as follows: in the next section we provide a reminder 
of the details of the obscrvables and their dependence on soft emissions as well as write 
down the resummed formula we derived in our previous work. In the following section 
we consider the general full NNLL result with variations of factorisation, renormalisation 
and resummation scales so as to derive the uncertainty on the perturbative result. Next 
we provide our comparisons to the experimental data and comment on the quality of the 
agreement as well as the potential need to include non-perturbative effects before providing 
a concluding discussion identifying future developments. We also provide for convenience 
an appendix where we list the main formulae we use in this work. 

2. The (p* variable and its resummation 

In this section we remind the reader of the main features of the resummed result computed 
in Ref. [28]. There we derived the dependence of the ot and (j)* variables on soft and, 
optionally, collinear gluon emissions from the incoming partons. In both cases one found 
that, in the soft limit, the dependence on emissions was essentially via a single component 
of gluon transverse momentum which was the one normal to the axis defined by the nearly 
back-to-back leptons. Here, and for the rest of this paper, we shall focus on (p*, which is a 
measure of the deviation of the Icpton opening angle A(p in the transverse plane from its 
value at Born level, Acp = ir. The (p* variable is defined as 



where (pacop = tt — A(p is the acoplanarity angle which vanishes at Born level, sin 6* derives 
from the angle in a boosted frame such that the leptons make angles 6* and n — O* with the 
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beam , M is the mass of the lepton pair and (pi is the angle of the gluon i with respect to 
the lepton axis in the transverse plane, with kTi the magnitude of its transverse momentum 
with respect to the emitting (incoming) partons. Requiring (p* to be equal to some fixed 
value thus involves the constraint 

6 (^<P* - " ^/°° dbMcos{bM4>*)'[[e'''^y^ , (2.2) 

where on the RHS we note that the sum over gluon emissions appears in factorised form 
in impact parameter (6) space as a product over emissions. When combined with the 
factorisation of the matrix element squared in the soft-collinear limit one can exponentiate 
the single-emission contribution (with account of running coupling) in b space. This results 
in the resummed form for the (p* distribution which reads [28] 

^((p*,M,cosd*,y) = ^ f dbM cos (bM^* ) e'^^^'^^) x S (xi , , cos ^* , 6, M) , 
d(p* sNc Jo 

where 

xi2 = ^e^y and b = ^. (2.3) 
VS 2 

The expression above is yet to be integrated over the dilepton invariant mass M, the scat- 
tering angle 9* and rapidity of the dilepton system (or equivalently the Z boson rapidity) 
y. In our previous work we took into account the experimental cuts over M and 9* but 
we integrated our expression over the full rapidity range. In this paper we will present the 
results in different rapidity bins, to able to compare to the data. 

The function R{b) in the exponent in Eq. ( p. 3D represents the single gluon emission 
contribution taking into account the running of a^- To be precise one can write 

R{b) = Lg^'^ (asL) + g^^^ (asL) + ^5^^) + . . . (2.4) 

TT 

where L = In (PM^). The functions Lg^^\ g^"^^ and ^g^^^ are the leading, next-to-leading 
and next-to-next-to-leading logarithmic contributions respectively. Further sub-leading 
terms are not shown, as they are beyond our current accuracy. The quantity S is basically 
the Born level result with the modification that the argument of the parton distribution 
functions (pdfs) entering therein is set as 1/6^ rather than a scale of the order of M^, the 
dilepton invariant mass squared. This modification of pdfs is due to DGLAP resummation 
of logarithms arising from hard collinear emission from the incoming legs and coincides 
with the standard treatment for the Qt variable. Also included in S are leading-order 
(O (as)) coefficient functions that are convolved into the pdfs. The precise form of T, can 
be found in Ref. [28] and in Appendix |A| for convenience. The radiator R{b) turns out to 
be precisely the same as for the Qt variable [28] and hence the only difference from the Qt 
case is the presence of the cosine function here rather than the standard Bessel function. 

In the current paper we present numerical results including full NNLL accuracy for 
R{b), whereas in our previous work we had kept only the leading term in as of the NNLL 

^Our definition of (j>* corresponds to the variable 0^ in the D0 study. 
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function g^^^ ~ a^L. The reason we have now used the full form of g^^^ is that, due to 
recent work [9] for the Qt variable, the complete result for g^^^ has become available, 
whereas previous results had assumed (in the absence of an appropriate calculation) that 
one of the coefficients involved in g^^\ which pertains to terms starting at order a^, would 
be the same as for the resummation of threshold logarithms. This turns out not to be the 
case and hence we use the recently computed result, which however does not have a visible 
numerical impact on our final results. 

Finally, an important issue that we handle in the next section is that thus far we have 
set all ambiguous scales in the problem to be of the order of the hard scale of the process, 
M. In reality our result is arbitrary beyond the NNLL accuracy of our resummation and 
the NLO accuracy of our fixed-order result in the standard way for any truncated fixed- 
order estimate. To assess this uncertainty is of course important before addressing the 
data and drawing conclusions, for instance, on the role of higher-order and sub-leading 
logarithmic terms omitted from our treatment and certainly before making any conclusive 
statement on non-perturbative effects. In the following section we thus also consider the 
role of scale variations on our resummed results for factorisation and renormalisation scales 
as well as varying the argument of the logarithm by a factor in a way that one generates 
sub-leading terms beyond NNLL, whose size can therefore be naturally estimated. 

3. Numerical evaluation, fixed-order matching and results 

The function R{b) can be found in Appendix ^ where we list the results for g^^\ g^"^^ 
and g^^^ as a function of the single logarithmic variable A = /3oas(Af ) In (PM^) . They 
are meaningful only in the range < A < 1, with the upper limit corresponding to the 
position of a divergence associated to the Landau pole, where at the one-loop level b = 
Aqcd^^- The lower limit at which R{b) vanishes is instead outside the jurisdiction of 
resummation. To handle these difficulties we need to limit the range of the b integral (see 
for instance Ref. [20]), such that 6min < b < 6max with bmm = 2M^^e^^^ and 6max = 
2M~-'^e~'''-^ exp (l/(2/3oas)). For our our current work, in the region below b = bmin we 
simply set the radiator to zero while for b > 6max one has no contribution since we assume 
the radiator to be infinite here. The precise choice we make for 6max corresponds to a cut-off 
on perturbative dynamics in the vicinity of the Landau pole and we note that in practice 
increasing 6max beyond 1/(3Aqcd)) the value we adopt here, has a negligible impact on 
our result. 

In principle of course this is not a rigorously justified procedure and one is free to 
choose other prescriptions to regulate the divergence of the perturbative result and replace 
it with genuine non-perturbative effects. In practice one may view the dependence on the 
precise position of 6max as a sign of the sensitivity to non-perturbative dynamics and in 
our particular case varying the position of 6max in the vicinity of the Landau pole does not 
change our numerical results significantly. Also related to this issue is the evolution of the 
pdfs down to the scale 1/6. In principle, we should evolve the pdfs down to l/bmax- In 
practice, we freeze the pdfs at a scale of the order 1 GeV, which is the minimum value for 
which the pdfs of the set we consider are tabulated. However, by evolving the pdfs with the 
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code HOPPET [29], we have checked that pushing the freezing scale below this value has 
negligible effects. We shall comment more on non-perturbative effects in the concluding 
section of this article. 

To obtain results which can be compared to the data over a wide range of (j)* , one 
needs to combine the resummed results with those from fixed-order codes which compute 
the (f)* distribution up to order a1, i.e. at the NLO level. In this context we note that our 
NNLL resummed result guarantees control over all logarithms that can arise up to NLO 
accuracy including the least singular a^L term. Having full control over all the logarithms 
at the two-loop order allows for a very simple matching formula 

da \ - ('^^\ + f ^ i'^^^ (3 1) 

'^'^* / matched V / resummed V ^<^* / fixed order V d*/"* / expanded 

where one simply combines the resummation with the fixed-order and subtracts the ex- 
pansion of the resummation to order to remove any double counting. We note that for 
the matching to be considered successful at small (j)* the expansion of the resummation 
and the fixed-order result should cancel and the matched result should thus follow the 
pure resummed curve while at large (p* the resummation should largely cancel against its 
expansion (up to relatively small order terms) and the result should tend to the pure 
fixed-order curve. 

Fig. |l] shows our basic resummed result with full NNLL resummation matched to NLO 
predictions from MCFM for the (ff distribution normalised to have unit area^ along with a 
comparison to pure fixed-order predictions from MCFM, with the (p* range being that over 
which the D0 collaboration has collected data [25]. The curves have been obtained for pp 
collisions at 1.96 TeV with the CTEQ6m set of parton densities [30], with the value of the 
strong coupling taken from the fit, as{Mz) = 0.1179. We consider leptons with invariant 
mass 70 GeV< M < 110 GeV, transverse momenta Iti,2 > 15 GeV and pseudorapidities 
|??i^2| < 2 in the central rapidity bin of the vector boson |y| < 1. This corresponds to the 
set of cuts adopted by D0 for muons. It is noticeable that the matched result and the 
fixed-order predictions agree at large (p* while the matched result behaves exactly like the 
resummation at smaller (j)* and the pure fixed-order result is seen to grow significantly at 
low (p* reflecting the logarithmic divergences contained therein. We have checked that the 
inclusion of full NNLL accuracy in R(b) has a negligible impact on our previous predictions 
for the same quantity, where we included only the first term of the NNLL function g^^^ [28] . 

We note that a difference from the Qt distribution is the absence of a peak in our final 
resummed result. As we have observed before this is related to the fact that two distinct 
mechanisms can produce a low (p* value: inhibition of gluon radiation which results in a 
Sudakov form factor and the vectorial cancellation between emissions that contribute to 
the one-dimensional sum in Eq. ( |2.2D . The latter mechanism which corresponds to the 
small-6 region of the b integral dominates the Sudakov behaviour at small values of p* so 
as to wash out the Sudakov peak in the (p* distribution which simply rises to a constant 
value. 



^We have checked that this normahsation differs from the NLO inclusive rate obtained with MCFM by 
O (a^) terms, as one would expect. 
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Figure 1: Figure illustrating the full NNLL resummed result to NLO for the (ff distribution. Also 
shown for comparison is the fixed-order result from MCFM while the range of 0* chosen is that 
over which data has been collected [25]. 



Having presented and discussed our basic resummed result and its matching to fixed- 
order we now turn to the question of assessing the uncertainty on the matched resummed 
result, as manifested by a dependence on the various arbitrary perturbative scales con- 
tained in the full result, which in the discussion above had all been set equal to the pair 
invariant mass M. The study of perturbative uncertainty is an important step as, while 
one should observe a significantly smaller scale uncertainty at low (jj* for NNLL resumma- 
tion as opposed to NLL resummation, one would nevertheless expect it to be important in 
the eventual comparisons to data and would certainly wish to take it into account before 
reaching conclusions on the size of non-perturbative effects. 

3.1 Perturbative uncertainties 

Here we shall deal with the issue of the perturbative uncertainty afflicting our calculation. 
In general, one can consider the strong coupling to be evaluated, as usual, at some renor- 
malisation scale which, although of the order of M, is not the same. Very similarly, the 
factorisation scale /x^r, which sets the scale where the pdfs are evolved from, ought not to 
be exactly equal to M. Another arbitrary scale in our perturbative calculation enters the 
argument of the logarithms we are resumming, ln(6^M^) — )• ln(6^//Q), where we refer to 
//Q as the resummation scale. For analogous studies for the Qt spectrum see, for instance, 
Ref. [21]. Keeping the full dependence on those scales the resummed expression in Eq. ( |2.3| ) 
becomes 

((/>*, M, cos = ^ r dbM cos{bM(P*)e-^'''^'^^'^Q'^^^ 

d(p* sNc Jo 

xT,{xi,X2,cos9*,b,M,ij,Q,fiR,fj.F) ■ (3.2) 
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Figure 2: Study of the perturbative uncertainties in case of NLL+LO (on the top left) and 
NNLL+NLO (on the top right). The black lines represent the variation of the resummation scale 
M/2 < /ig < 2M, while the red ones the variation of renormalisation and factorisation scales 
Af/2 < {^jf ^fip} < 2M, with 1/2 < ^h/M-F < 2. The orange (yellow) band is obtained by varying 
all scales independently, requiring the ratio of any two of them to be between 1/2 and 2. At the 
bottom the uncertainty bands are normalised to the curve with all scales set to the pair mass 



Explicit formulae are reported in Appendices ^ and The resummation is then matched 
to a fixed-order calculation, which also depends on renormalisation and factorisation scales. 
The dependence of the resummed and matched result on these arbitrary scales is one order 
higher than the accuracy we are working at, i.e. it affects terms which are at least N^LL 
and NNLO. Thus, varying them around the dilepton invariant mass M provides us with 
an estimate of the size of those perturbative contributions which are beyond our accuracy. 
In doing so our theoretical prediction becomes a band, and the central value, where all the 
scales are set equal to each other and to the pair mass, does not have any special physical 
meaning. 

In Fig. Q we perform a study of the theoretical uncertainty that affects our result. As 
an example we consider the set of cuts for the muons in the central rapidity bin \y\ < 1, 
but similar results can be found by looking at the other bins. We study two different 
levels of accuracy: NLL+LO, on the left, and NNLL+NLO, on the right. The black lines 
represent the variation of the resummation scale M/2 < < 2M, while the red ones 
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Figure 3: Comparison of the theoretical prediction NNLL+NLO for the 4>* distribution to the 
experimental data collected by the D0 collaboration in the case of muons, in two different vector 
boson rapidity bins, \y\ <1, on the left and 1 < jyl < 2 on the right. 



represent the variation of renormalisation and factorisation scales M/2 < < 2M, 

with 1/2 < fJ-n/fJ-F < 2. The orange (yellow) band is obtained by varying all scales 
independently, requiring the ratio of any two of them to be between 1/2 and 2. In the same 
figure, at the bottom, the uncertainty bands are normalised to the curve with all scales 
set to the pair mass. We note that at low (p* the dominant source of uncertainty is given 
by the fig variation, i.e. it comes from sub-leading logarithmic terms in the resummation. 
We also note that the uncertainty in the small (p* region is almost halved in going from 
NLL+LO to NNLL+NLO. In fact we have that the size of the band is O (20%) in the 
former case, while it is O (10%) in the latter. This is consistent with the NNLL+NLO 
uncertainty band found in the case of the Qt spectrum [21]. 

Another source of theoretical uncertainty comes from the parton distribution functions. 
However, we find that these effects mostly cancel once the distribution has been normalised 
to the inclusive rate, leaving an uncertainty at the percent level, much smaller than the 
band obtained with scale variations. 

Having estimated the theoretical uncertainty of our calculation, we can now compare 
it to the experimental data. This is the topic of the next section. 

4. Comparisons to data 

The D0 experiment at the Fermilab Tevatron performed a measurement of the (j)* distri- 
bution in bins of \y\, using 7.3 fb^^ of data [25]. Different kinematical cuts are applied 
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Figure 4: Comparison of tlie theoretical prediction NNLL+NLO for the ip* distribution to the 
experimental data collected by the D0 collaboration in the case of electrons, in three different 
vector boson rapidity bins, \y\ <1, on the top left, 1 < |y| < 2 on the top right and \y\ > 2 at the 
bottom. 



for muons and electrons. More specifically, the invariant mass of the leptons must lie in 
the range 70 GeV< M < 110 GeV. Moreover, muons must satisfy /ti,2 > 15 GeV and 
\rji^2\ < 2, while electrons must satisfy ^7^1,2 > 20 GeV and |?/i,2| < 1-1 or 1.5 < |r/i^2| < 3. 
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Figure 5: Comparison of our NNLL+NLO tlieoretical prediction for the (f)* distribution to the one 
obtained with the program RESBOS, in the case of electrons in the central rapidity bin \y\ < 1, on 
the left, and the forward one \y\ > 2, on the right. 

We now compare our theoretical calculations to the data, taking fully into account these 
experimental cuts. The theoretical predictions are normalised to the area under the curve, 
as are the data. In Fig. ^ we compare the data to our theoretical prediction for the muon 
channel, in the central rapidity bin |y| < 1 (on the left) and in the outer one 1 < |y| < 2 
(on the right). The yellow bands represent the theoretical uncertainty obtained by vary- 
ing the three scales ij,q, hr and jip independently, as explained in Sec. |3.1| . The solid 
(green) curve is the theoretical prediction obtained with all scales set equal to the pair 
mass = = fip = M. The plots in Fig. ^ instead show the comparison to the elec- 
tron data. In this case we have three different rapidity bins: \y\ < 1 (on the top left), 
1 < |y| < 2 (on the top right) and \y\ > 2 (at the bottom). The last rapidity bin for the 
electrons provides an opportunity to probe fairly small-x values (x < 10^^) with conven- 
tional resummation techniques. In this study we are mostly interested in the low cp* region 
of the distribution, hence we do not show the large (p* tail of the outer rapidity bin because 
the experimental errors, as well as the theoretical uncertainties, are rather large and the 
comparison is not very instructive. 

We observe, on the whole, within the scale uncertainties, an excellent agreement over 
a large range of (p* for both muons and electrons in all rapidity bins. This includes the 
lowest (p* values where one may expect non-perturbative effects to play some role. 

It is worth noting that we obtain a comparable description of D0 data to that provided 
by the commonly used event generator RESBOS [17, 18], however, with different under- 
lying physical assumptions. This is shown in Fig. ^, where we compare our NNLL+NLO 
theoretical predictions, in case of electrons in the central (left) and forward (right) rapid- 
ity bins, to the ones obtained with the RESBOS, run with the same event selection as in 
Ref. [25]. We note that the height of the plateau predicted by RESBOS is lower than ours, 
because of the presence of non-perturbative effects, modelled by the BLNY form factor'^. 



^The BLNY form factor [16] is primarily controlled by one parameter g2. We use the default value 
32 = 0.68 GeV^. 



There is a striking difference between the sizes of the theoretical uncertainties of the two 
predictions. Fohowing [25], the RESBOS uncertainty is obtained by varying fiR and 
simultaneously by a factor of two, although we do not include here PDFs uncertainties. 
This procedure does not capture the main source of uncertainty, which, in our approach, 
is given by variations of the resummation scale ^q. Finally, in the forward rapidity bin, 
we also we show the RESBOS prediction including small- a; broadening effects which were 
obtained by fitting RESBOS to semi-inclusive-deep-inelastic-scattering data [26]. The re- 
sulting distribution is broader and, consequently, its height is lower. Small-x broadening 
is disfavoured by the data, as discussed in [25]. 

Our analysis of the theoretical uncertainty has led us to the conclusion that it is hard 
to make any statements on the size of non-perturbative effects, except to set a range of 
values based on the upper and lower edges of our uncertainty band. Nevertheless, it is 
possible to take our central values, with all scales set equal to the pair mass, and correct 



it to the data by adding a non-perturbative smearing to the radiator Eq. (2.4): 



i?NP {bM) = R (bM) + 5Np6^ , (4.1) 

which corresponds to assigning a Gaussian-smeared intrinsic transverse momentum to the 
incoming partons. We can now look for the NP parameter (7np which gives the best 
description of the data. The results are shown in Fig. |6| and |^ with different choices of gt^p, 
inspired by the literature on Qt resummation. We note that the spread of results obtained 
by varying the NP parameter from g^p = GeV^ to g^p = 1 GeV^ is similar to the band 
describing the perturbative uncertainty. The best value of the NP parameter is different in 
different bins. We obtain (^np — 0.5 GeV^ for central rapidities and (7np — 0.3 GeV^ in the 
more forward regions. However, we stress that this procedure is misleading since we have 
no reason to assign any special role to the curve where all the perturbative scales are set 
equal to the dilepton mass. In doing so we would be ascribing the perturbative uncertainty 
to a universal NP parameter and the use of this parameter for studies of related variables 
could result in erroneous conclusions. 



5. Comments and conclusions 

In the current paper we have carried out a complete phenomenological study for the (p* 
variable recently measured by the D0 collaboration [25]. We have included resummation 
to the next-to-next-to-leading logarithmic level and matched the resummation to the full 
NLO calculation from MCFM [11]. We carried out a study estimating the theoretical 
uncertainty on the matched resummed prediction and provided comparisons to data. We 
also investigated the potential role of non-perturbative effects in the comparison of theory 
to experiment. 

One of the aims of our study was to exploit the fact that the measurement errors on (p* 
are significantly smaller than for the standard Qt variable [23], to study as precisely as pos- 
sible the role and accuracy of resummed calculations and eventually also non-perturbative 
effects. In this context, we noted that the pure fixed-order (NLO) estimate was inadequate 
to describe the (p* distribution over a very wide range of measured (p* values, indicating 
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fiV. |y|<1 



nV, i<lyl<2 




Figure 6: Distributions obtained with different values of the non-perturbative parameter 5np, in 
the different rapidity bins for muons. All curves are obtained setting fig = = fip = M. 



both the need for resummation as well an an opportunity to compare different resummation 
approaches (including parton shower event generators) with a view to identifying significant 
physics differences and further developing resummation tools and event generator models. 

We observed that our basic matched resummed results depend as usual on various 
arbitrary scales, specifically, the renormalisation scale fiji, the factorisation scale fip and 
the resummation scale . Varying these scales gives an idea of the impact of sub-leading 
terms omitted from our treatment and, as a result, we derived an uncertainty band for our 
prediction, which reflects the theoretical error on our estimate. 

We found that, within the aforementioned error band, we get an excellent description of 
the (j)* distribution with hardly any variation in the quality of our description in the different 
bins in y, the lepton pair rapidity. We would thus deduce that the size of non-perturbative 
effects cannot be ascertained unless one reduces further the theoretical uncertainty with, 
for instance, an N^LL resummation. We also directly studied non-perturbative effects and 
found that reasonable values of a non-perturbative Gaussian smearing generated a spread 
of results comparable in size to our perturbative uncertainty. We stress that taking our 
central value and correcting it to the data with a non-perturbative effect, while possible, 
is misleading since it leads potentially to ascribing physics, which could be of perturbative 
origin, to universal non-perturbative effects. A preferable approach to non-perturbative 
effects from our viewpoint would be to theoretically study them for the (p* variable using 
for instance the techniques based on a universal infrared finite extension of ag as applied 
for instance to Qt resummation [19] and the related case of energy-energy correlations in 
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eV, |y|<1 



eV, 1<|y|<2 




Figure 7: Distributions obtained with different values of the non-perturbative parameter gNP, in 
the different rapidity bins for electrons. All curves are obtained setting hq = fij^ — ij,p — M. 



e~^e annihilation [20]. We postpone this investigation to future work. 

As far as further developments are concerned, one could envisage, given the excellent 
agreement of the precise perturbative results with the (j)* data, that with a suitable fit 
range one may use the results presented here to fit for the strong coupling, ag, from the 
D0 data. On the more theoretical side one useful development would be to try and for- 
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mulate the resummation of the (p* variable in a 6-space-independent framework. The main 
reason for this is that the use of 6-space requires us to provide arbitrary prescriptions to 
deal with regions of integration which are not strictly under the control of resummation, as 
well as serves to potentially contaminate our results with non-perturbative effects via the 
use of parameters such as 6jnax5 designed to avoid the Landau pole. While we have noted 
the stability of our results against reasonable variations of such parameters, it is clearly 
desirable to remove them altogether, if possible. Resummed studies exist for the Qt distri- 
bution which do not employ the 6-space formalism, both within traditional resummation 
methods [32] and within the techniques of soft-collinear effective theory [9, 10] and it would 
be interesting to see if such methods could be used more generally for variables such as (p* 
and other similar studies. We note that similar methods to those used in our current study 
have also been used to study azimuthal decorrelations between coloured final state particles 
(specifically for jets see Ref. [33]), which represents an extension of the Qt resummation 
formalism to processes with colour in the final state. 

Lastly, since it is straightforward to extend our studies to the LHC, we would also 
advocate measurement of the (p* variable at the LHC and suggest that studying 0* together 
with the Qt distribution [34] would provide invaluable information especially when detailed 
comparisons are made to the variety of theoretical calculations and Monte Carlo models 
that are currently available, including our forthcoming results. 
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A. Resummed result 

The resummed expression for the cp* distribution has been reported in Eq. (2.3), with all 



the perturbative scales set equal to the dilepton mass and in Eq. (p.2|), where they were 
kept separate. The function S is given by 



"s(Mj?.) gCorr ( M 



T,{xi,X2,cos6*,b,M,fj.Q,fiR,fiF) = e " v''Q/ x 

T,{xi,X2,cos6*,b,M,fiQ,fiR,fiF) , (A.l) 



9 ( — ) = ^ — - -Cf In — . (A.2) 



where 

_corr / iH \ ^ h:£ 

This contribution ensures full control of the logarithms at NNLL in the expansion, that 
is the dependence on the resummation scale /ig starts at N'^LL. Note that S acquires a 
dependence on the impact parameter b because of the resummation of logarithms of b via 
DGLAP evolution, which then determines the scale of the parton distribution functions 



-15- 



embedded in S. Including the contributions from the Z as well as from the virtual photon, 
we have 

S = (1 + cos2 r) {Q\ - 2QqViVqXi + {A^i + ^')(^^ + Vq)X2) 

+ cos0*(-4QgAAxi + ^AiViAqVqX'i)T- , (A.3) 

where a sum over quark flavours q is implied. The above equation is naturally written as 
the sum of the two terms with different angular dependence: the first one is proportional to 
(1+cos^ 0*) and represents the parity conserving piece of the electro- weak interaction, while 
the term involving cos Q* is the parity violating piece. We notice that upon integration over 
the full Q* range, as well as over symmetric intervals, the parity violating term vanishes. 
The coefficients Ai^q and Vi^q are the electroweak couplings for lepton / and parton g, 
explicitly given by: 

Af = Tf and Vf = r| - 2Q/ sin^ Bw, f = l,q, (A.4) 

where Tj is the third component of the isospin. We also have introduced 

M2 - M|) 
" ''(M2-M|)2 + r|M|' 

2 

_ V2GfMI 



Ana 

In Eq. (A.3), are explicitly given by 



(A.5) 



= (C ^ h)qixi,b, HQ, flR, flp) {C i2)q{x2,b, flQ, fiR, flp) 

± {C0 fi)g (xi, 6, i^Q, fiR, fip) (C (g) f2)g (2:2, 6, fJ-Q,f^R, fJ-p)- (A.6) 

The convolutions involving the matrix of coefficient functions C and the vector of 
parton densities fi^2 for incoming hadrons 1 and 2 respectively can be explicitly written as 

(C fOq f^Q,f^R, f^F) = ^Cqo, («. (^) ^ ^ ' ' ^^'^^ 

where i = 1, 2 and a sum over all flavours a is implied. The combinations of scales in the 
argument of the strong coupling and of the pdfs are such that, when b = bmm = /^g^ and, 
consequently R = 0, they reduce to hr and hf, respectively. 

The coefficient functions Cqa represent perturbative corrections to the collinear branch- 
ing of an incoming parton a to a parton q which annihilates with q to form the Z boson. 
We note that the collinear enhanced terms generated by such a branching are incorporated 
to our accuracy into the scale of the pdfs fi via their dependence on the impact parameter 
b. Thus the coefficient functions represent only the non-logarithmic constant terms. Be- 
cause of their collinear origin, the scale of the coupling constant in the coefficient function 
is naturally of the same order as the one for the scale which pdfs are evolved to. However, 
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there is in principle a hard (process-dependent) contribution as weh, coming from virtual 
corrections, characterised by kx ~ M. We consistently choose to work in the so-called 
"Drell-Yan scheme" in which one sets this contribution to one [8] . 

The coefficient function admits the following expansion in the strong coupling constant: 



In this work, we consistently resum next-to-next-to-leading logarithms in the radiator. 
When the exponential is expanded, the O (a^) contribution to coefficient function generates 
terms which are of the same order as the ones in g^^^ and hence, in principle, should also be 

included. However, these contributions start at O (a^) and hence wc would expect their 
effect to be numerically as significant as the one due to the coefficient A^"^^ , which wc have 
found to be small. Therefore, in the current paper wc use the first order approximation of 
the coefficient functions and we leave the inclusion of the O (a^) contributions to future 
work. The explicit form of the O (a^) coefficient functions is 

C« (x, ^ =Cf{^^- 4) (5(1 - x) + (1 - x) + In , (A.9) 



where pf) are the LO DGLAP splitting functions. 



B. The radiator to NNLL 

In this appendix we collect the explicit expression for the functions gi, which enter into the 
radiator: 

+%(3)fo.L,^,^), (B.l) 



where L = ln(6^//Q) and = agifJ-R) is the MS coupling. The explicit expressions for the 
functions g^"^^ are 
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A(^^ A + In (1 - A) 

7r/3o A 

ln(l - A) + 



(B.2) 



/3)(A) 



vr^o 



A 



7r2/32 \^i_A 
A + ln(l-A) 1,2 



+ ln(l - A) 



1 - A 



+ -ln^(l-A) 



AW 



A^ + lna-A)")!!!^, 



A2 A(2)/3i / A(3A - 2) (1 - 2A) ln(l - A) 



(B.3) 



27r2^2 (1 _ ;r/3o 1 - A ttP^ \ 2(1 - A) 



{1-xy 



(Pf 1-2A u,Ml-A) 



1- A 



(/3o/32(1-A) + /3i2a) 



+ 



A 



2(1- A)2 
bW(3i ( A 



(^0/92(2 -3A)+;3f A) 



1-A 



+ ln(l-A)l+4^^1n'-| 



2fe 



1-A 7r/3o (1 - A) 



/ A 1-2A , ,^ 



In 



2 ' 

(B.4) 



with A = OLs{liR) PqL. 

The coefficients A^*) and B^*) are given by 



(B.5) 
(B.6) 



.rs^ C'f /^2 /^245 134 2 11 4 22 \ 1 / 418 40 ^ 56 \ ^, 
^(3) = (72 ^2 _^ 4 _^ ^ _|_ ^ Cs (B.7) 

16 V V 6 27 45 3 V 2 ^ V 27 27 3 ^7 ^ ' 



+ ^Cpuf ( -rr + 16C3 ) - ^n} ) + lnPod2 , 



55 

Y 



2 /^Q 



(B.8) 



■77 vr^ + ^ 



+ A(^)ln^. 



(B.9) 



The coefficient ^2 has been recently determined [9]: 



d2 = CfCa ( ^ - 28(3 



112 

Yn 



CpUf . 



(B.IO) 
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We have also introduced the coefficients of the QCD /3-function 
llCA-2n/ 



Po = 



127r 

-»2 



247r2 



2857/ 



^2 = 64;^3 ■ 
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